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Q_|| Abstract. We consider the Cauchy problem for the nonlinear Schrodinger equa- 

in space dimensions d > 3, with initial data close to 
^n , a sum of A'^ decoupled solitons. Under some suitable assumptions on the spectral 

On I structure of the one soliton linearizations we prove that for large time the asymp- 

totics of the solution is given by a sum of solitons with slightly modified parameters 
and a small dispersive term. 

0. Introduction 

O 

g ; In this paper we consider the nonlinear Schrodinger equation 

2; (0.1) iipt = -A^ + F{\^\'^)^, {t,x) eRxR'^, d>3. 

Q^i For suitable F it possesses important solutions of special form - solitary waves (or, 

shortly, solitons): 

e^'^ifix -b{t),E), 



^ I where w, 7 G M, f , c G M'^ are constants and (/? is a ground state that is a smooth 

■ positive spherically symmetric, exponentially decreasing solution of the equation 

(0.2) ~Aip + Eip + F{ip'^)ip = 0. 

Solitary wave solutions are of special importance not only because they are simple 
and sometimes explicit solutions of evolution equations, but also because of the 
distinguished role they appear to play in the solution of the initial value problem. 
This is best known for completely integrable equations like the cubic Schrodinger 
equation 

(0.3) i^t = -i^xx - IV'I'V'- 

In the general position case the solution of the Cauchy problem for this equation 
with rapidly decreasing smooth initial data has in L2(M) the asymptotic behavior 

N 



d ■ 1 

' ^ = ujt + ^ + -X ■ V, b{t) = vt + c, E = uj + ^ > 0, 
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where Iq = —d^, and /+ is some function in L2{^). The number A*", the function 
/+ and the soHton parameters {jj,Ej,Vj,Cj) depend on the initial data. Due to 
possibihty of explicitly integrating equation (0.3) with the help of inverse scattering 
methods, they can be described by effective formulas in terms of initial condition. 
See, for example, [22] for these results. 

Numerical experiments have shown that even without the presence of an inverse- 
scattering theory, solutions, in general, eventually resolve themselves into an ap- 
proximate superposition of weakly interacting solitary waves and decaying disper- 
sive waves (see [11], for example). While exact theory confirming the special role 
of solitary waves as a nonlinear basis with respect to which it is natural to view 
the solutions in the limit of large time is not generally available, partial indication 
is provided by stability theory of such waves. A considerable literature has been 
devoted to the problem of orbital stability of solitons following the work of Ben- 
jamin [1], see also [7, 13, 14, 29, 34, 35]. The problem arises in connection with the 
Cauchy problem for equation (0.1) with initial data of the form 

(0.4) V|t=o = ^(^'^o) + XO' 

where xo is small in the Sobolev space H^{R'^). It was shown that under certain 
additional conditions the solution 'tp{x,t), t > remains close (again in the space 
H^iR'^)) to the surface 

This notion of stability establishes that the shape of the wave is stable, but does 
not fully resolve the question of what the asymptotic behavior of the system is. 

The first asymptotic stability results were obtained by Soffer and Weinstein in 
the context of the equation 

(0.5) #t = -A^ + [y(x) + A|^r-i]V', 

(see [27, 28] and [30, 31, 32, 33, 36] for the further developments related to this 
model). The solitons for (0.5) arise as a perturbation of the eigenfunction of the 
operator —A -|- V{x) and, in contrast to the case of equation (0.1), they have a 
fixed center, which simplifies the analysis to some extent. 
For the one- dimensional equation 

(0.6) i^t = -Vx^ + F{\^\^)^ 

the asymptotic stability of solitons was studied in the works of Buslaev and author 
[4, 5]. We considered the Cauchy problem (0.6), (0.4) and proved that in the case 
where the spectrum of the linearization of equation (0.6) at the initial soliton has the 
simplest possible structure in some natural sense, the solution tp has an asymptotic 
behavior of the form 

= e'^+^{x-b+it),E+) + e-''°'f+ + o{l), ^+=uj+t + ^+ + ^, b+=v+t + c+ 

as t ^ -\-oo, where the parameters {'j^, E^,V-^C-^) of the limit soliton are close 
to the initial ones (0,i?o,0,0) and /+ is small. Some asymptotic results in the 
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framework of significantly freer conditions on the linearization were obtained in [5] , 
see also [6]. Recently the analysis of [4, 5, 6] was extended to the multidimensional 
case (0.1) by Cuccagna [8, 9]. 

As a natural generalization of the above situation one can consider the case 
of several weakly interacting solitons. Assume that one has a set of solitons 
^ijJaj+t 2 (^(,x — Jjq-^ Eqj), j = 1, . . . , N, that are well separated either in the orig- 
inal space or in Fourier space: for j ^ k, either or niin|6"^(t)| is sufficiently 

large, where Vjj^ = voj — vok, b^ki'^) — ^oj — ^ofc + '^jk^. In the second case we shall 
assume that the "collision time" t'^-u = — ''^^^0^,2''° is "bounded" from above, see 

subsection 1.4, (1.7) for the exact formulation. 

Consider the Cauchy problem for equation (0.1) with initial data close to a sum 

N 

i=i 

If all the linearizations constructed independently from the solitons (p{EQj) satisfy 
the spectral conditions introduced in the case of one soliton, one can expect that 
as t ^ +00 the solution V' looks like a sum of N soliton with slightly modified 
parameters plus a small dispersive term. In [23] this was proved in the case d = 
1, N = 2, see also [19] for the asymptotic stability results for the sums of solitons 
in the context of KdV type equations. The goal of the present paper is to extend 
the result of [23] to the multidimensional case d > 3 (omitting also the restriction 
N=2). The main new ingredient in the analysis is a combination of the estimates 
for the linear one soliton evolution obtained by Cuccagna in [8] with the ideas of 
Hagedorn [15]. 

The structure of this paper is briefly as follows. It consists of two sections. In 
the first section we introduce some preliminary objects and state the main result. 
The second contains the complete proofs of the indicated results, some technical 
details being removed to the appendices. 

1. Background and statement of the results 
1.1. Assumptions on F. Consider the nonlinear Schrodinger equation 
(1.1) iVt = -AV' + i^dV'HV', it,x) eRxR'^, d>3. 

We assume the following. 

Hypothesis HO. F is a smooth function, F{0) = 0, F satisfies the estimates 

nO>-Ce^ \F^"HO\<Ce-", a = 0,1,2, 
where C > 0, ^ > 1, q < ^, p < 

Set g{0 = E^ + F{e)C 

Hypothesis HI. 

(i) There exists ^0 > such that g{^) > for ^ < Co; 5(0 < C > Co (ind 

g'iCo) < 0. 

(ii) There exists ^1 > such that J^^ dsg{s) = 0. 

Further assumptions are given in terms of the function 

/(C,A) = -A^<?'(0 + (A + 2)5(0. 
We consider ^0 of (HI) and assume: 
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Hypothesis H2. For any C > Co there exists a A(C) > 0, continuously depending 
on I, such that I{t, X) > for < t < ^ and I{t, X)<0 fort> ^. 

We suppose hypotheses (HI, 2) to be true for E in some open interval A C M+. 

Under these assumptions equation (0.2) for E A, has a unique positive spheri- 
cally symmetric smooth exponentially decreasing solution ip{x, E), see [2, 20]. More 
precisely, as |x| — >■ oo 

^{x,E)^Ce-^\-\\x\-^. 

This asymptotic estimate can be differentiated any number of times with respect 
to X and E. 

We shall call the functions w{x, a) = cxp{i(3 + iv-x/2)(p{x~b, E), a = (/3, E, b, v) 
G ]^2d+2 soliton states. w{x,a{t)) is a solitary wave solution iff a{t) satisfies the 
system: 

(1.2) p' = E-^-^, E' = 0, b' = v, v' = 0. 

1.2. One soliton linearization. Consider the linearization of equation (1.1) on 
a soliton w{x,a{t)): 

iXt = (-A + Fi\w\^))x + F'{\w\^)i\wW + w^x)- 
Introducing the function /: 

/=(j), x(x,t)=exp(z$)/(y,t), 

^ = m + ^, y = x-b{t), 

one gets 

ift = L{E)f, L(E) = Lo(E) + ViE), L^{E) = (-A + E)a3, 

V{E) = Vr{E)a^ + iV2{E)a2, Vx = F{^'') + F'(^2)^^ V2{E) = F'(^2)^^ 
Here (72, CT3 are the standard Pauli matrices 

-i\ (I 



''^ = \i J' "'-yo -1 

We consider L as an operator in L2(M'^ ^ C^) defined on the domain where Lq is 
self adjoint. L satisfies the relations 

CTsLas = L*, aiLai = —L, 

where "^i ~ q^- "^^^ continuous spectrum of L{E) fills up two semi-axes 

{—oo,E] and [£',00). In addition L{E) may have finite and finite dimensional 
point spectrum on the real and imaginary axis. 
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Zero is always a point of the discrete spectrum. One can indicate d + 1 eigen- 
functions 

io = '^(_\) ' ^3 = 'fiyj (J) ' i = 1' ■ ■ ■ d, 
and d + 1 generalized eigenfunctions 

L^j = 0, L^d+i+j =^j, j = 0,...,d. 

Let M be the generalized null space of the operator L. Under assumptions (HO, 1,2), 
the vectors ^j, j = 0, . . . ,2d + I, span the subspace M iff 

see [34, 20, 8]. 

We shall assume that 

Hypothesis H3. The set Ao of E e A such that 

(i) zero is the only eigenvalue of the operator L{E), and the dimension of the cor- 
responding generalized null space is equal to 2d + 2; 

(ii) ±E is not a resonance for L[E); 
is nonempty. 

Obviously, the set is open. 

Remark. zizE is said to be a resonance of L{E) if there is a solution ip of the 
equation {L{E) =F E)ip = such that < x >~* ^ G L2 for any s > 1/2 but not for 
s = 0. zLE can never be a resonance if d > 5, see lemma A4.3. 

Consider the evolution operator e~**^. One has the following proposition. 

Proposition 1.1. For E e Aq and any xo,xi € M'^, 

(1.3) II {x - xo)-"" e-^^(^)*P(£;)/||2 < C {t)-"'^ II {x - x,r fh, ^0 > ^, 

where P{E) is the spectral projection onto the subspace of the continuous spectrum 
ofL{E): 

KerP = M, RanP = (asM)-^. 

The constant C here is uniform with respect to xo,«i € M.^ and E in compact 
subsets of Aq. 

This proposition is an immediate consequence of the Lp- Lg estimates of e~*^*P 
proved by Cuccagna [8]. For the sake of completeness we sketch the proof of (1.3) 
in appendix 4. 

1.3. The nonlinear equation. We formulate here the necessary facts about the 
Cauchy problem for equation (1.1) with initial data in H^{R'^). 
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Proposition 1.2. Suppose that F satisfies (HO). Then the Cauchy problem for 

equation (1.1) with initial data ip{x,0) = ipo{x), ^ H'^iJR.'^) has a unique solution 
ip in the space C(]R — >■ H^), and ip satisfies the conservation laws 

where C/(^) = dsF{s). Furthermore, for allteR 

WmWm <c(||V'o||//i)||V'o||//i, 

where c : M+ — »• M+ is a smooth function. 

The assertion stated here can be found in [10, 11], for example. 

1.4. Description of the problem. Consider the Cauchy problem for equation 
(1.1) with initial data 

N 

(1-4) V'|t=o = V'o e if^ nLi, tl^o = ^w{;aoj) + Xo, 



(1.5) aoj = {l3oj,Eoj,boj,voj), min \vjk\ >vo>0. 

Here f^^. = voj — vok- Set bf^,^ = boj - bok, j k. Write 6^^ as the sum 

hO . .,0 

l^-Oj "jk - ^jk ^jk^jk^ ^jk ^jk — S'fc ~ L,0 12 ■ 

\Vjk\ 

For j k we define the effective small parameter ejk- 

■ (min|60,(t)| + ^0,1)-!, if t% < k < r% >, 

jfcl 



(1-7) e,k 



\v^Au\ otherwise. 



where 6°fc(t) = 6°^ + tv^jk^ k is a fixed positive constant. 
Assume that 
(Tl) e = maxejk is sufficiently small^; 

(T2) Eo, €Ao,j = l,...,N. 

Our goal is to describe the asymptotic behavior of the solution ip as t ^ +oo, 

provided xo is sufficiently small in the following sense: 
(T3) for some m' , ^ + ^ = 1, m > 2p + 2, | + 2<m< + 2 if d > 4 , 

4 < m < +2 if d = 3, the norm 

•A/" = llxolli + llxollm' 

is sufficiently small. 

Here xo stands for the Fourier transform of xo • 
Our main result is given by the following theorem. 



"Sufficiently small (large)" assumes constants that depend only on wq , ft and Eoj ,j = 1, . . . ,N. 
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Theorem 1.1. Fort > the solution o/(l.l), (1.4) admits the representation 

N 

where \Ej{t) - EQj\, \vj{t) - VQj\, j = 1,...,N, ||x(t) ULani™ are small uniformly 
w.r.t. t>0, and as t ^ +oc, 

||x(t)IU = o(t-"^*-^)). 

Moreover, there exist vectors a+j = {P-\-j,E^j,b+j,V-^-j), such that as t ^ +oo, 

\a,{t)-a+,{t)\ = Oit-'), 
for some S > 0. Here (J+j{t) is the trajectory of {1.2) with the initial data (T+j(0) = 

2. Proof of the theorem 

Up to some technical modifications the main line of the proof repeats that of 
[23]. 

2.1. Splitting of the motions. Following [23] we decompose the solution if) as 
follows. 

N 

(2.1) i^{x,t)=Y,w{x,aj{t)) + x{x,t). 

i=i 

Here (Jj{t) = {Pj{t), Ej{t),bj{t),Vj{t)) is an arbitrary trajectory in the set of ad- 
missible values of parameters, it is not a solution of (1.2) in general. 

We fix the decomposition (2.1) by imposing the orthogonality conditions 

(2.2) (^fj{t),asik{Ej{t))^=0, j = l,...,N, k = 0, . . . ,2d + I. 
Here 

fj = ( j'.) ' X{x, t) = exp{i^j)fj{yj,t), 

$j = Pj{t) + Vj ■ x/2, yj = X- bj{t), 

< •, • > is the inner product in L2(M'^ ^ C^). 

Geometrically these conditions mean that for each t the vector fj (t) belongs to 
the subspace of the continuous spectrum of the operator L{Ej(t)). 

For ^ of the form (1.4) with min(|w°j,.| + |6°jt.|) sufficiently large, and with xo 

sufficiently small in some Lp norm, the solvability of (2.2) is guaranteed by the 
non-degeneration of the corresponding Jacobi matrix, see lemma Al.l. So, one can 
assume that the initial decomposition (1.4) obeys (2.2). To prove the existence of 
a decomposition (2.1), (2.2) for t > 0, one can invoke a standard continuity type 
argument, see appendix 1 for the details. 



Rewriting (2.1) as an equation for % one gets 
(2.3) ixt=H{a{t))x + N, 

Here 



N 



V{w) = {F{\w\') + F'{\w\')\w\')as + F'{\w\')(^_l, ) , w,=w{x,aj). 
The nonlinearity N is given by the following expression 

N 

N = TVo + 5^e-3*^/(a,-)fo(2/i,£^,-), 



No = F{\i^, + x\ 



N . . . . N 



'J ' 2 

where 7j, cj are defined as follows. 

/?,■(*) = dsiEjis)-^-^^- \ ' )+ljit), h^{t) = dsv^{s) + c^{t). 

In terms of parameters (7, E, c, v) (1.2) takes the form 

y = 0, E' = 0, c' = 0, v' = 0. 

Substituting the expression for xt from (2.3) into the derivative of the orthogo- 
nality conditions, one gets for j = I,... ,N 

ieiEj)E'. = {Nj,ase'^^ioi- - bj,Ej)) + (^f^,l{aj)io{Ej)) , 



(2.4) eiEj^ = {Nj,ase'^^ id+ii- - bj,E,)) + {fj,liaj)i,+,{Ej)) , 



in{Ej)c'^ = - ((iV,-,(73e^*^e«+i+fe(- - hj,Ej)) + (^fj,l{aj)id+i+k{Ej] 



k=l,...,d 



9 



Here 



=^0+ Yl ^K)X + Yl e'"^*'^K)fo(yfc,i^fc), j = i,...,N, 



d 2 1 II ii2 

e = ^ll<p||2, n=2ll^ll2- 



The right hand side of (2.4) also contain the derivative a', which enters hnearly 
in /(cTfc). In principle, system (2.4) can be solved with respect to derivative and 
together with equation (2.3) constitutes a complete system for a and x- 



(2.5) 



iXt = H{d{t))x + N{d,x), 



(2.6) 



a' = G(a, x), x\t=o = Xo, <7j(0) = aoj. 



2.2. Integral representations for x- In this subsection we follow closely the 
constructions of Hagedorn [15] (developed in order to prove the asymptotic com- 
pleteness for the charge transfer model), see also [21]. We start by rewriting (2.5) 
as an integral equation 



(2.7) xit)=Uoit,0)xo-i fuoit,s) 

Jo 



N 



Y^V,{s)xis) + N 



ds, 



Here Uo{t,T) = e»(*-^)^<^3^ y. = V{wj). 

Next we introduce the one soliton adiabatic propagators Uf{t,T): 

iUft{t,T) = Lj{t)Uf{t,T), Uf{t,T)\t=r = I, 

Ljit) = -Aa3 + Vj{t) + Rjit), Rjit) = iro,(t)[Pj(t), P,(i)]To^(t), 
V,(i) = Tojmm{Eoj)T*{t)T*.{t), P^it) = Tjit)P{Eoj)T*{t). 



Here 



4 / JO 



ds{Vj{s) -Voj), 



(Joj{t) = {/3oj{t), EQj,boj{t),voj) being the solution of (1.2) with initial data (Toj (0) 
aoj- Obviously, 



where 



Pfit) = Toj{t)Pjit)T^jit). 



Write the solution x as the sum: 
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Using the adiabatic evolution Uj^{t, r) one can write the following representation 
for hj{t) 

{2.8) hj{t)=Uf{t,0)Pf{0)xo-i f Uf{t,s)Pf{s)[ V Vm{s)x{s) + D^{s)]ds, 

Jo , ■ 

Here 

(2.9) Dj=N + {Vj-Vj)x-RjX- 
Combining (2.7), (2.8) one gets finally 

(2.10) X = (I) + (!!) + (Ill) + (IV), 

where ^ 

(I)=Uo{t,0)xo-iYl [ dsUoit,s)Vj{s)Uf{s,0)Pf{0)xo, 
. Jo 

(II) = - J] f dsK,{t,s)Vm{s)m, 

i.m. Jo 

(III) = -i f dsUo{t,s)D, 
Jo 

(IV) = -5^ f dsK^{t,s)D^{s). 

Here 

(2.11) D = N + Y, (V,% + (V,- - V,)x) , 
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Kj{t,s) = ^ dpUo{t,p)V,{p)Uf{p,s)Pf{s). 

The relations (2.4), (2.7), (2.10) make up the final form of the equation which is 

used to prove theorem 1.1. 

2.3. Estimates of solitons parameters. Following [4, 23] we consider (2.4), 
(2.7), (2.10) on some finite interval [0, ti] and then study the limit ti +oo. On 
the interval [0,ti] we introduce a natural system of norms for the components of 
the solution if): 

N 

Moit) = J2 ITj (*) - 0oj\ + \Ej{t) - Ejo\ + \cj{t) - boj\ + \vjit) - voj\, 

M,it) = Y,\\<y,>-''xm2, M2(t) = ||x(t)||2p+2, ^>^, 
without loss of generality one can assume that m = 2p + 2. 
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These norms generate the system of majorants 

Mo(t) = sup Mo(t), Mi{t) = sup Mi{t)p-''^{t), 1 = 1,2, = Mk{h). 

0<T<t 0<T<t 

Here 1</Lti<|if(i = 3 and 1 < //i = ^ for d > A, 112 = — 2p^2 ) ' 



tjk being "the coUision times" that are defined as follows. We set tjk = if t^f^ < 0. 
For (j, k) such that i"^ > 0, we define tjk by the relation, 



Jo \v%\' 



^0 

where 

I Vjk(ti), \it> ti, 

Let us mention that 

(i) tjk are well defined provided \vjk{t) — Vj^.\ < vq, < t < ti; 

(ii) the conision times tjk belonging to the interval [0,ti] "do not depend on ti\ 
It follows directly from the definition of Mq that 

(2.12) \e'j{t)l \a'j{t)\ < Mo{t)+Mi{t), \bj{t) - bj{t)\ < Mo{t), 

(2.13) \<I>j{x,t) -^j{x,t)\ < Mo(i) <x-bj{t) > +Mo{t) [ ds\c'j{s)\, 

Jo 

where 

bj (t) = boj (t) + aj (t) , ^j ix,t) = (3oj (t) + Oj (t) + voj ■ x/2. 
It is also easy to check that bjk = bj — bk admits the estimates 

(2.14) \bjkit)\>c\v%\\t-tjk\, 

(2.15) \bjkit)\ > c(min \b%is)\ + \v%\\t - tjk\) - c, t% <k< r% > 

provided Mo{t) < c for < t < ti. Here and below c is used as a general notation 
of positive constants that depend only on Vq, n and eventually on Ej,j = 1, . . . , N, 
in that case they can be chosen uniformly with respect to Ej in some finite vicinity 
of Eoj. 

Consider relations (2.4). Since 

INK II in-ui n^/ W' + W'"^' ifi>>i^ 

\No\<c\2_^\wj\\wk\[l + \x\) + < I ,2 .f / 1 

I \x\ if ^' < 2 



j,k lAI " ^ 2 
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and 

j ^ k, Xj = {'yj,Ej,cfj,Vj), bjk = bj — bk, one gets immediately from (2.4) 

(2.16) \Xj{t)\ < W^(M)[J]e-'=l''"'WI + {ml{t) + Ml{t)) p''^'{t)]. 

i,l 

We use PF(M) as a general notation for functions of Mq, Mi, M2, which are 
bounded in some finite vicinity of the point = 0, Z = 0, 1,2, and may acquire 
+00 out some larger vicinity. They depend only on vq, kq, Ejo, j = 1, . . . ,N and 
can be chosen to be spherically symmetric and monotone. In all the formulas where 
W appear it would not be hard to replace them by some explicit expressions but 
such expressions are useless for our aims. 

Combining (2.13), (2.16) one gets 

(2.17) \^j{x,t) - ^j{x,t)\ < W{M)Mo{t) <x- bj{t) > . 
Integrating (2.16) and taking into account (2.14), (2.15) we obtain 

(2.18) Mo < W(M)[e + M?+Mi]. 

Consider the vectors kj{t) = {I - Pf{t))x{t), kj{x,t) = E?=o^ kji{t)e'^^'^^^i{x- 
bj{t), Eqj). The orthogonality conditions (2.2) together with (2.12), (2.17) lead 
immediately to the estimate: 

(2.19) \kji{t)\ < M/(M)Mo(t)||e-'=l"-^^ WIx(*)l|2 < W {M)Moit)Mi{t) p"' (t) . 

2.4. Linear estimates. To study the behavior of solutions of the integral equa- 
tion (2.10) we need some estimates of the evolution operators U^{t,T)P^{T). The 
necessary estimates are collected in this subsection, the proofs being removed to 
the appendices. 

Lemma 2.1. For any xq, xi G R'^, < t < t < ti, 

(2.20) II {x - xo)-''°Uf{t,T)Pf{T)fh < WiM) {t - t)-"!^ \\ {x - xip fh- 

The function W here is independent of xq,xi and ti. 

See appendix 2 for the proof. 
Remark. Due to the representation 

Ufit,T)PfiT)f = Pfit)Uoit,T)f-ij\sUf{t,s)Pf{s){Vj{s) + Rjis))^^^^ 

and the estimate 

\{R,{t)f){x)\ < I^(M)e-^l=^-''^«l(|^;.| + |a;.|)||e-'=l--''^(*)l/||2 
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(2.21) < W^(M)e-^l^-^^(*)lMo(i)||e"'=l''"''^'(*^l/||2, 

(2.20) leads immediately to the inequality 



(2.22) II {x - xo)-"" Uf{t, T)Pf{T)f\\2 < W{M)- 



Pl + WfWpO 



where 2 < pi < < P2 < oo, — + ^ = 1, i = 1, 2. Obviously, the same estimate 
is valid for Kj{t,T): 



(2.23) II {x - xo)-"" K,it,T)f\\2 < W{M)- 



p{ + lU IlPo 



The key point of our analysis is the following lemma that is essentially lemma 
3.6 of [15]. 

Lemma 2.2. Introduce the operators Tjki{t, t), j,k,i = 1, . . . ,N, k 

Tjki{t,T) = Aj{t)Kk{t,T)Ai{T), 
where Aj{t) is the multiplication by < x — bj{t) . Then, for <t <ti 



fdr\\TMt,T)\\ < W{M){e^l+Mo{t)), 
Jo 



with some ui > 0. The norm \\ ■ \\ here stands for the L2 L2 operator norm. 

See appendix 3 for the proof. 

2.5. Estimates of the nonlinear terms. Here we derive the necessary estimates 
of D, Dj. We write D as the sum: 

D = D^ + D^+D^, 

where 

j 

A^oo = F(|V'.P)(_^4) -E^(K-|')(_"^.) +V(^.)X-EV,-X' 



X 



In a similar way, 

Dj = D''+D^ + D', j = l,...N, 
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where 

k 

Estimating A^oo by 

(2.24) |iVoo|<c(l + |x|)5;K-|K|, 

and using (2.12), (2.17), (2.19), (2.21) one gets 

< PF(M)[(1 + Ixl) J]e-'=(l^-'''l+l^-'"=l) 

+ ^e-c|x-6,|(|^.|^Mo(t)|x|+Mo(OMi(i))]. 

i 

which together with (2.16) leads to the inequahty 

(2.25) \\D^L,nL„ P°||LxnL. < W(M)[e-=l''^'=WI + (MqMi + + Mi)p'^^ (t)]. 
Consider D^, D^. We estimate them as follows. 

\D^ + D^\ < W^(M)[|V^,||xP + IxP + IxP^+a if = 3, 

(2.26) \D'\ < ^^(M)|V^,||xP, \D^\ < \xf^+\ if ^ < p < 1, 

+D'^\<WiM)\x\'^P+\ if_p<l. 
These inequalities imply for r' = 

\\D^ + D^h^nL^, < M^(M)[M2 +mJ"^M| +M2^^]p''i(t), if d = 3, 

(2.27) 

\\D^L,nL^, + \\D^L^,nL^, < W(M)[M? + Mr^M| +M^+>^Hi), if \ <P<'^ 

Ili^' + ^'llvnL^, < W(M)M^+VHi), if ^- 

2.6. Estimates of x in -C'2,/oc- To estimate Mi{t) we use representation (2.10). 
By (2.22), for the first term (I) one has 

(2.28) II < (i)||2 < wim)j\f < t >-'^/^ . 

Consider expression (II): 

(2.29) II < yj (II)||2 < W{M)M^{t)Y^ [ ds\\Tjki{t, s)\\p^^ (s). 
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By lemma 2.1, 

\\Tjki{t,s)\\<W{M)<t-s>-''/^ . 
So, the integral in the right hand side of (2.29) can be estimated as follows. 



^ ds\\Tjkiit,s)\\p^^{s)< {^j^ ds\\Tjkiit,s)\\p'"\s)^ ' (^j^ ds\\Tjki{t, 



s) 



2f»l 



< T^(M)(M^ ^ ^^-^K/ ds<t-s >-'^/2 //2(5)^ 

< 6 = 1 — z/2 = Ovi- At the second step here we have used lemma 2.2. Thus, 

(2.30) II < y, >-^ (II)||2 < VF(M)(M^ + er|)Mi(t)p^Hi)- 

Consider the two last terms in the r.h.s. of (2.10). By (2.25), (2.23) (with 

Pi = 2, p2 = oo) one has 

II < yj >-'Uo{t,s)D°{s)\\2, II < yj Km{t,s)Dl{s)\\2 < W{m) <t-s>-'"^ 

(2.31) x[J]e-l^^'=(")' + (Mo(t)Mi(t)+M2(t) + M2(t))p^i(s)]. 

Using (2.27), (2.23) one can estimate the contribution of D^, D'^ as follows. 
\\<yj >-''Uo{t,s){D\s)+D\s))h, || < y,- K„,{t,s){D\s) + D\s))h 

(2.32) <W{M)^l{t)+Ml'{t)]\t- s\-^'^ <t-s>-'''+''^ p^'{s). 

Here 1 < ri = 1 + min{p,p~^} < 2. Combining (2.31), (2.32) and integrating with 
respect to s one gets 



< 



yj (III)||2, II < % >-^ (IV)||2 < W{M)[Y, [ ds 

i Jo 



i^k 



<t-S >'^/2 



+(MoMi +M^^)p''i(i)], 

or taking into account (2.14), (2.15), 
(2.33) 

II < y, (III)||2,|| < y.,- (IV)||2 < VF(M)[€ + MoMi + M2 + M^^]p^^(t). 

Combining (2.28), (2.30), (2.33), one obtains 

Ml < W{M) [M + e''^ + M^Mi + + M^^ ] . 
Changing if necessary the coefficient function W one can simplify this inequality: 
(2.34) Ml < W{M)[N + e''^ +M2^\. 
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2.7. Closing of the estimates. Here we derive a estimate of x which will 

close the system of the inequalities for the majorants. To estimate L„ - norm of % 
we use representation (2.7). By (2.24), (2.27), 

||Ar||„, < W(M)[5]]e-^l^*'=WI +M?+M^^)p^i(i)]. 

k,i 

As a consequence, 

(2.35) M2 < iy(M)[A/' + e^-^2 ^Mi]. 

Here we have made use of the inequality 



/ ds- ; — < W(M) 

Jo ~ 



<t-tik >^'^ 



which is an immediate consequence of (2.14), (2.15). 
Combining (2.18), (2.34), (2.35) one gets 

(2.36) Ml, M2 <W{M){M + Mo <W{M){M^ + 6^"^), 

1/3 = min{i, 1/2,1 — /i2} > 0, the coefficient functions PF(M) being independent 
of ti. These inequalities mean that for TV and e sufficiently small M can belong 
either to a small neighborhood of zero or to some domain whose distance from zero 
is bounded from below uniformly with respect to Af, e. Since are continuous 
functions of ti and for ti = are small only the first possibility can be realized. 
This means that for J\f and e in some finite vicinity of zero, 

Mi(t), M2(i) < c(A/' + e'^3), Mo{t) < c{M^ + e^''^), Q<t<ti. 

The constant c here is independent of A/", e, ti. Since ti is arbitrary these estimates 
are valid, in fact, for alH > 0. More precisely, one has 

(2.37) 

Mo(i)<c(Ar2 + e2-3), M,{t)<c{M + e''')p^{t), M^it) < c{M + e''')ff^{t), 
where Poo(i) is the weight function corresponding to ti = 00: 

j.k 
j^k 



*Tk Vjk{s) ■ v% _ 

"'^ 1-0 12 ~ jfe' 



i°°fc = 0iftO, <0,and 
if t% > 0. 

By (2.15), (2.16), the estimates (2.36) imply the existence of the limit trajectories 
{t) = W+3 it) , E+j ,b+j{t),v+j), j = 1,...,N, 

b+j {t)=v+jt + b+j , v+j = voj + / dsvj (s) , 

^0 
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|2 



ds{c'j{s) + Vj{s)-v+j), 
/3+,{t) = iE+, - + E+j = Eoj + dsE'^is), 

OO 



Obviously, as t ^ +oo, 

|i;,(t)-i;+,| = 0(t-2^^+i), 
\bj{t) - b+j{t)l \Pj{t) - P+j{t)\ = 
Appendix 1 

Here we outline the arguments needed for the proof of the existence of a decom- 
position (2.1) satisfying (2.2) for all t>0. We begin with the following lemma. 
Given N solitons w{aQj), a^j = {Poj, EQj,boj,voj), j = 1,...,N, we define the 
effective coupling parameter 5{ao), uq = (ctqi, . . . , (Joat), 

For X e Lp{R'^), a = (ai, . . . , ajv), aj = {Pj,Ej,bj,Vj) € M x ^ x M'^ x R-^, 
j = 1,. . . ,N, consider the functionals Fj^i{a, x', ^o), j = 1; • • • ) -^j / = 0, . . . 2d + 1, 



^ - \ 

+ ^w{aok) -w{(Jk),a3Ci{aj) j , 

k=i I 



where 



w 



{^^, Q{x,a) = e'^'^-+'^'^^ii{x-b,E), a = iP,E,b,v) 



Set Fj — {Fj^o, . . . , Fj^2d+i), F — (Fi, . . . , Fn). 

Lemma Al.l. Let E^j G Aq, j = I, . . . , N. There exist constants uq > ^, 5q> 0, 
K > depending only on Eqj, j = 1, . . . ,N such that ifS{ao) < Sq and \\x\\p < '^o 
then the equation 

F{d,x;^o) = 

has a unique solution a = (at, . . . ,(Tn), a being a function of x, that satisfies 



1 



(Al.l) \pj - boj + -{vj - voj) ■ boj\ + \Ej - Eoj\ + \bj - boj\ + \vj - voj\ < K\\x\\ 



v 



Remark. It follows directly from (Al.l) that 
(i) for some constant Ki 



N 

\\x + ^w{aQk) - w{ak)\\p < Ki\\x\\p, 
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(ii) if for some pair (j, k), t^-^ < < f% > then the new colHsion time tjk = —^^7-^ 
satisfies a similar estimate with a constant k = k;o(1 + 0(||x||p)- 
Proof of Lemma Al.l. Let us pass from a to a new system of parameters A = 
(Ai, . . . , Aat), 

= (03 - Poj + ^{vj - voj) ■ boj,Ej,bj - boj,Vj - voj). 

We represent F{a, x', <?o) as the sum 

F = F'^ + F^ + F^, 
F^ = ^Xj,Eoj),^ = {^o,...'^2d+i), 

^1= E (CoKfc)-CoK),d(c^,))- 

k, k^j 

At last, 

= GiiXj, fj), x{x) = e'^^^+'^^fjix - boj), 

Gi{X,f) =< f,asCi{^) > is a function of / and A. 
The direct calculations give 

(A1.2) |detVA$(A,£;)||,^(3_^3_„^ = e\E)n''\E). 

Set Ao = (Aoi,...,Ao7v), \oj = (0,£oi,0,0). By (Al.2), 

AT 

(A1.3) |detVxF°||x=Xo = \{e\Eo^)n^\Eoi) 

is nonzero if E^j E Aq, j = I, . . . ,N . 

Consider F^. It is not difficult to check that for A in some finite vicinity of Aq 
the derivative V^-F^ satisfies the inequality 

(A1.4) IV^F^I <C5(ao), 

constant C depending only on E^j. 

By the implicit function theorem, the desired result is a direct consequence of 
(A1.3), (A1.4). ■ 

To prove the existence of a decomposition (2.1) satisfying (2.2) for all t > we 
use some standard continuity type arguments. Since V' G C(M ^ if^) there exists 
a small interval [0, ii] where the constructions of lemma Al.l can be used. This 
leads to a representation (2.1) satisfying the orthogonality conditions for t G [0, ti\. 
For the components of such a representation estimates (2.15), (2.36) give 

\Eaj -E\< C{M^ + e''^'), {\vjk\ + \bjk\)-^ < Ce, 
\\xmm<C{J\f + e''^), 

which allows us to extend decomposition (2.1), (2.2) on a larger interval [0, ti + 12] 

with some ^2 > 0. On this new interval the same estimates hold, so one can 
continue the procedure with steps of the same length t2- As a result, one gets a 
decomposition (2.1) satisfying (2.2) for all t>0. 
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Appendix 2 
Here we prove lemma 2.1. Consider the equation 

(A2.1) ixt = C{t)x, C{t) = {-A + E)a^+V{t) + i[P'{t),P{t)l 

V{t) = T{t)V{E)T*{t), P{t) = T{t)P{E)T*{t), 

where T(i) = -Be(t),a(t),o- We denote the corresponding propagator by U{t,T) 
Clearly, 

We shall assume that for some positive constants n, R, 6i, 
(A2.2) \e\t)\ + \a' {t)\ < n, 



(A2.3) \e"{t)\ + \a"it)\ < ^ < - ti) > 



1=0 



t G R+. Here LsN, = io<^i<---<iL- One has the following lemma. 
Lemma A2.1. For any Xq, Xi G , t > Q,t > 0, 

II {x - bo)-''° U{t,T)PiT)fh <C{t- r)-"/' II {x - x,r fh, 

provided n is sufficiently small and R is sufficiently large: n + R~^ < C. 

In this appendix we use C as a general notation for constants that depend only 
on M, 6, E and can be chosen uniformly with respect to E in compact subsets of 

Aq. 

It follows from (2.12), (2.14), (2.15), (2.16) that for M in some finite vicinity of 

zero the functions 6j, aj satisfy assumptions (A2.2), (A2.3) with 5i = 2/ii — 2, ti, 
I = 1,. . . ,L, being the collision times tik, i,k = 1, . . . N, i ^ k. So, lemma A2.1 
implies lemma 2.1. 

Proof of lemma A2.1. Lemma A2.1 follows from proposition 1.1 by a simple 
perturbation argument. On the intervals [t;, / = 0, . . . L — 1 we introduce the 
following linear approximations 0\t), a\t) of 9{t), a{t): 

ft 



e\t) = e{t) - j^^ dsj^^ ds, (^1 - o"{s,) 

fti+i rti+i _ J. 

- / ds dsMr^ — tW"{si), 

a\t) = a{t) - j'^ ds £ ds, - r,{^^^)^ a"{s,) 



rti+i rti+i _ 

/ ds i dsir]{- -)a"{si). 

Jt Js H+i - ti 
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1 forlCl < i, 
forl^l > |. 



Here r] G C°°(M), r/(0 = 

For t & [tL, 00) we define the corresponding 9^~^^{t), a^+^(i) as follows. 



ds y dsie"{si), 
ds J dsia"{si). 



Clearly, for t G I = 0, . . . ,L, tL+i = 00, one has 

(A2.4) \e{t)-e\t)\, \a{t)-a\t)\<CR-'^, 

On the interval [t;, t^+i] one can pick out the leading term of (A2.1) in the form 
(A2.6) ixt = C\t)x. C\t) = (-A + E)a^ + V\t), 

V\t) = T'{t)V{E')T'\ T\t) = i?AHt),aHt),r^ 

A'(t) = e\t) = -r, E =E + -r- A 



2 ' ' 4 

We denote the propagator corresponding to (A2.6) by U\t,T). Clearly, 

u\t,T) = r'(t)e-^(*-^)-^(^')r'*(T), p\t)u\t,T) = u\t,T)p\T), 

where P'(t) = T\t)P{E^)T^' (t). 

Consider the expression x{t) = U{t,T)P{T)f, ti < t < ti^i. 

For ti < t < tij^i we write as the sum x = /i + A;, = P\t)x{t). Since 
x{t) = P{t)x{t), the 2d + 2 dimensional component k is controlled by h: 

||e^l^-«WI A;(i) II2 < C{\e{t) - e\t) I + \a{t) - a\t) | + |r' | + 1^ - |) ||e-^l^-"«lx(t) II2 



(A2.7) < C{R-^ + n)||e-^l"-"«lx(t)||2 

for some 7 > 0, provided n, are sufficiently small. In the last inequality we 
used (A2.4), (A2.5). 

For h one can write the following integral representation 

(A2.8) h{t) = P\t)ho{t)-i dsP\t)U\t,s)[V\s)ho{s) + R\s)x{s)l 
where 

R\t) = Vit)-V\t) + i[P'{t),Pit)]. 
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Obviously, 

\V{x,t) - V\x,t)\ < C\e{t) - e^{t)\ + \a{t) - a\t)\ + |r^| + \E - £;'|)e-^l^-"WI 



(A2.9) < C{R-^ + n)e-^l^-"WI, 



(A2.10) \[P'{t),P{t)]f\ < (7ne-^l"-"(*)l||e-^'"-"WI/||2. 

Estimates (A2.7), (A2.9), (A2.10) and representation (A2.8) together with propo- 
sition 1.1 imply immediately that ioT ti < t < t < and for any ^ G M the 
following inequality holds 

II {x-a{t))-''°x{t)h<t-T + ^>''/' 
(A2.11) < C sup f II {x - a{s))~''° e'^^"^^'-^') P{T)f\\2 <s-t + ^ y"^/^) , 

T<S<t ^ ' 

where C do not depend on ^. (A2.11) implies in particular, that 

(A2.12) II (x - a(t))-'^" C/(t, r)P(r)/||2 <C<t-T >-'l'' \\ {x - xoY" /lb, 

XiyeW^,ti<T <t< ti+i, l = 0,...,L. 

To prove that this estimate is in fact true for any < r < t we use the induction 
arguments. Assume that one has (A2.12) for r < t < < oo. We need to show 
that then the same is true ior r < ti < t < t^+i. For t G we write 

U{t,T)P{T)f = U{t,ti)U{ti,T)P{T)f. Using (A2.12) and the representation 

U{t,T)P{T)f = e'^t-r){^-E)a,p^^y 



(A2.13) -i y* dse^(*-")(^-^)'^^ (V(s) + P(s)]) U{s, T)P{T)f, 

one checks easily that 

II < X - a{t) >-'° e^^'^^(*-*')i7(i;, r)P(T)/||2 <C<t-T y-"/^ \\ {x - xo)"" fh- 

By (A2.11), this implies that (A2.12) is valid for < r < t < and thus, for 
any < t < t. Moreover, by (A2.13) one can replace a{t) in the left hand side of 
(A2.12) by any Xi G R'^: 

\\{x-Xl)-''°U{t,T)P{T)f\\2<C\\{x-Xor fh, 0<T<t. m 
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Appendix 3 

Here we prove lemma 2.2. We start by proving a similar result for the "free" 
operators Tj*,,^ , t) : 

T%i{t,T)=Aj(t) dpUo{t,p)Vk{p)Uoip,s)MT). 
Lemma A3.1. For i ^ k, t > 0, one has 

(A3.1) fdr\\T^,S,r)\\<WiM)e'f^ 

Jo 



with some ui > 0. 
Proof. Since 

one has 



\TUt,T)\\<C<t-T> 



-d/2 



(A3.2) I{t)= f dT\\T%S,r)\\<C 

Jo 



< t > 



In this appendix the constants C depend only on Eok. 

For t > 2/9, where p is a small positive number, we write the integral I{t) as a 
sum of two terms I{t) = Io{t) + Ii{t), 



m = dT\\T^^^{t,T)i 

TP^,{t, r) = Aj{t) f ' dsUoit, s)Vk{s)Uo{s, t)Mt), 

Jt+p 



Ii{t) being the rest. Obviously, 

(A3.3) h{t) < Cp. 

Consider Jq {t) . To estimate this expression we write Tj*^^ [t, r) in the form 

/ n-ll n-12 \ 

rpp / jki jki \ 

■^jki — 1 _7'21 _'7-22 I > 

\ jki jki J 

where 



Jt+p 

Jt+p 



jki 



r+P 



(A3.4) rf,\{t, T)f = T^\{t, T)f, T^\{t, T)f = T^]{t, T)f 
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Consider Tjj^^{t,T). Since Hilbert-Schmidt norms dominate operator norms, we 
have 

\\T^\{t,T)f<C J dxdy<x>-''^<y>-^^ \B]ki{t,T)\\ 

where 



B]u{t,T)= [ dsit - s)-'/\s - r)- 

Jt+P 



-d/2 

It+p 

I i|3;-z + bj(t)-6t.(s) i|z-y + b^.(s)-b^(T)| 

X / dze Vi{z)e ^^^^ 



Integrating by parts in the second integral and taking into account (2.13) one gets 
immediately the estimate 

\B]ki{t,r)\ < Wim)i< x> + <y >)p-' 

p 

j ds{t - S^'^S - t)-'^/2 < ^.^.^^^ g^^) 

where 



t-p 

X 



djki{t,s,T) = ' + 



(t-s) {s-rY 

Here the function W do not depend on p. As a consequence, one has for < a < 
min{l,f -i}, \bjk{t)\ + \bik{r)\>0, 

t 

||^-fei(i,T)|| < W^(M)p-l-'^+2« j ds<t-S >-<^/2+«< s - T >-''/2+« 

T 

(A3.5) X<t-T >-'^/2+2a _ ^.^(^)| + _ ^)|^.^(i)|)-a. 

Here we made use of the obvious inequality 

J ds<s >-"< s-py-" \dis + dal"" <C <p {\di\ + \d2\)-", 

M 

provided a>l,0<a<l, (ii,d2Sl^^,C being independent of di, d2- 
Integrating (A3. 5) and taking into account (2.14,15) one gets finally, 



(A3.6) / ' dT\\T^\{t,T)\\ < W{M)p-'-''+^-e^. 

Jo 

In a similar way, one has for Tj^^-{t, r) 

\\rlk'{t,r)f <C I dxdy < X >-'^< y >-^^ \B%,{t,r)?, 
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.t-p 



Jt+p 

dze e2'*(^+'"=(^)'^)y2(2)e ^(^) 



which imphes 



t-2p * * 







\\T^ki{t,T)\\ < W{M)p-^-'^+^" J drJdsKt-s >-'«/2+« 

r 



(A3. 7) < W{M)p-^-'^+^"efk. 

Combining (A3.2), (A3.3), (A3.4), (A3.6), (A3. 7) one obtains 

Iit)<WiM)ip + p-^-^+^^efk), 

which leads immediately to (A3.1) with i/i < 2+d-2a ' ' 
Let us introduce the operators Tjj^-{t,T): 

T^„i{t,T) = A,{t) dsUo{t,s)%{s){I - P^{s))Uo{s,t)Mt). 

It is not difficult to check that for any a < 1, 

\\A,{t)Uo{t,s)%{s){I-P^{s))U,{s,T)MT)\\ 

< W{M) <t-s >-''/2< s-T >-'^/2+«< bikir) >-'^ . 

As a consequence, 
(A3.8) 

/ dT\\Tli{t,T)\\<W{M) f dT<t-T>-''/^+^<kkiT) >-«<W(M)ef^. 
Jo Jo 

At the last step here we have used (2.14), (2.15). 

Proof of lemma 2.2. This lemma follows directly from (A3.1), (A3. 8) and the 
following representation 

(A3.9) -i^ dp jysAjit)Uoit,p)Vk{p)P^ip)Uiip,s)Rk{s)Uo{s,T)Ai{T) 
(A3.10) -i jypAj{t)Uo{t,p)Vk{p)Pk{p)Ak\p)TLiP,r) 
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(A3.11) 

" dp r dsA,{t)Uo{t,p)%{p)P^{p)U^{p,s)[%{s) + Rk{s)]A^\s)T^,,^^^ 



T J T 



We estimate the right hand side of this representation term by term. Using lemma 
2.1 and inequahty (2.21) one gets 



(A3.12) / (iT||(^3.9)|| < VF(M)Mo(t). 

Jo 

Expression (A3. 10) can be estimated as follows 
(A3.13) 

/ dT||(^3.10)|| < W{m) [ dr [ dp<t-p ||rO,,(p,r)|| < W{M)ef,. 

Jo Jo Jt 

In a similar way, 

f dT\\{A3.11)\\<W{M) f dr f dp T ds < t - p >-'^/^ 

Jo Jo Jt Jt 

(A3.14) <p_5>-'^/2 \\T^,.{s,T)\\<W{M)el. 

Combining (A3.1), (A3.8), (A3.12), (A3.13), (A3.14) one gets lemma 2.2. ■ 

Appendix 4 

Here we discuss the proof of proposition 1.1. Since only the weighted estimates 
are needed, rather then follow [8, 37, 38] we use the approach of [16, 17, 18]. It turns 
out that the arguments of [16, 17, 18] can be applied almost without modifications. 
So, we describe only the main steps of the proof, referring the reader to [16, 17, 18] 
for most of the details. 

We start be recalling briefly some basic properties of the free resolvent Ro{X) = 
{-A + E-X)-'^ 



, . _ ,^ 1 . Let i/*'* stand for the weighted Sobolev 
—{-A + E + X)~^J 

spaces: 

H*'^ = {/, \\f\\H^,s ^ II < X >^ (1 - A)*/V||2 < oo}. 

We denote by B(H^'*, iJ«i'*i) the space of bounded operators from H^'* to 
Set L| = B{H''*) = BiH^'^H"'^). If s > 1 and t G M the resolvent Ro{X) 

which is originally defined as B{L2) valued analytic function of A G C \ ( — cjo, —E] U 
[E, oo) can be extended continuously to the C+ = {im A > 0} when considered 
as a S(iJ*'*, iJ"*'*"*"^) valued function. The following properties of Ro{X) are well 
known, see [16, 17, 18, 37, 38] and references therein. 

Lemma A4.1. Let A; = 0, 1, ... . If s > k + 1/2, then the derivative Rq'\x) G 
B{H^''^,H-^'°) is continuous in XeC+\{E, -E}, with 

(AAA) 4'=)(A) = 0(|A|-('=+^)/2)^ 

in this norm as X ^ oo in C+. 

The behavior of Ro{X) for A close to zizE is described by the following lemma, 
see again [16,17,18]. 
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Lemma A4.2. As X ^ E, Ro{X) admits the following asymptotic expansion in 
For m odd: 

I I 
(A4.2) iio(A) = G„o(A - Ey + G,,i(A - i^^+i + 0((A - EY^'), 

3=0 j=0 

for m even: 

I I 
(A4.3) RoiX) = 5^ G,- o(A - Ey + ln(A - E)^, Gj,i{X - Ey + o((A - E)'), 
j=o j=0 

where I = 0,1,..., s > C{l,d), the coefficients Gj^k belong to 

Gj^i = for j < if d is odd and for j < if d is even. Representations 

(A4.2), (A4.3) can he differentiated with respect to X any number of times. 

Here (A — and ln(A — E) are defined on the complex plane with the cut 

along [E, oo). The explicit expressions for the constants C(Z, d) can be found in [17, 
18, 19]. Similar expansions hold as A ^ —E. 

For A € [E, oo), consider the operator 

/ + i?o(A + iO)V : L^' L2 ^ 

s > 1. 

Lemma A4. 3. Let E ^ Aq. Then KeT{I + Ro{X + iO)V) is trivial. 

Proof. We start by the case X = E. Let ■0 £ Ker(/ + GqV). This implies that ip 
belongs to L2(M'^)+ < x >-('^-2) Loo(R'^) and satisfies 

Lij; = Etp. 

Hypothesis H3 then allows us to conclude that -0 = 0. 

We consider next the case X > E. Let ip G Ker(/ + Ro{X + iO)V). Since V is 
spherically symmetric, one can assume that ip{x) = f(r)Y(Lo), r = \x\, lo = 

f G L2{R+;r'^-^ < r >-^' dr) and Y G LalS'^"^), 

As.d-iy = //„y, /x„ = n(d-2 + n), 
for some ra G {0, 1, . . . }. Then / has to satisfy 



(A4.4) /„/ = 



(-—-'^—^— + E+—)u + V 
dr'^ r dr 



/ = A/, 



/'(O) = if n = 0, /(O) = if n > 0, 

and as r ^ 00, 

(A4.5) f = cr-^H';}\kr)(^\+0{e-'^^), 7 > 0, 
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for some constant c. Here k = (X — E)^^^ > 0, v = n+ Hi^'' is the first Hankel 
function. Asymptotic representation (A4.5) can be differentiated with respect to r 
any number of times. 
The Wronskian 

wif,9) = r'^~\< f',9 >R2 -<f,9' >R2) 

does not depend on r if / and g are solutions of (A4.4). Calculating w{f,f) one 
gets 

2ik\c\^ = 0, 

which implies that ^ G L2. Since E ^ Aq, this means that ^ = 0. ■ 
Consider the full resolvent R{X) = (L - X)'^. R{X)P {R{X)) is a ^(La) 
valued holomorphic (meromorphic with the only pole in zero) function of A G 
C \ (—00, —E] U [E, 00). R{X) satisfies the relations 

(A4.6) aiR{X)ai = -R{-X). 

The analytic properties of R{X) near the cuts {—oc,—E], [E,oo) are collected in 
the two following lemmas. In both of them wc assume that E G ^o- 

Lemma A4.4. For s > 1, R(X)P can be extended continuously to C+ as a 
B{L2, L2^) valued function. Moreover, if s > k + ^ then R^''\X)P exists and 
continuous for X G C+ \ {E, —E} and 



(A4.7) i?(^)(A)P = 0(|A|- 



(fe + l)/2N 



in 5(1,2, L2 *) as X ^ 00 in C+. 

Lemma A4.5. As X ^ E, R{X) admits the following asymptotic expansion in 
B{L2, L2'''). 



For m odd: 



l-i 



(A4.8) i?(A) = %o(A - Ey + %i(A - Ey+h + 0((A - E)% 

j=0 j=0 

for m even: 

I CO 

(A4.9) R{X) = YI1 Bj,k{^ - - + +«((^ - ^)')' 

j=0 fc=0 

where / = 0, 1, . . . , s > C(/, d), Bj^k G ^(-^i, L2'), S,,fe = /or /c = 1, j < 

if d is odd and for k > if d is even. Representations (A4.8), (A4.9) can be 

differentiated with respect to X any number of times. 

These results is a standard consequence of the corresponding properties of the 
free resolvent (lemmas A4.1,2) and lemma A4.3, see [16, 17, 18]. 

Consider the propagator e~**^. Lemma A4.4, together with (A4.6), (A4.8), 

(A4.9) allows us to represent the expression ^e~**^P/, 5^, f,g & Cq°{R'^) in the 
form 

(A4.10) (e-'^'^Pfg) = dAfe"^^* (^(A)/,5) - e'^* {S{X)a,f,a^g)], 
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where 

^(A) = ^{R{X + iO) - R(X - iO)). 

It follows from (A4.8), (A4.9) that as X ^ E, £{X) admits the following asymptotic 
expansion in B{L2,L2^) with s sufficiently large. 
For d odd: 

(A4.11) £{X) =£o{X- E)"^ +0{{X- E)^), 

for d even: 

^ r 0(ln(A -E){X- Ef) if d = 4, 



(A4.12) £{X) = £o{X - E)— + 



0{{X-E)i) ifd>6. 



£q G B{L2,L2^). These expansions can be differentiated with respect to A any 
number of times. 

Combining (A4.10), (A4.7), (A4.11), (A4.12) one gets immediately [18] 

II < X e-'*^Pf\\2 <C <t >-'^/2 II < a; /II2, 

provided s is sufficiently large. To recover proposition 1.1 it is sufficient now to 
inject this inequality in the following representation for e~**^P 

Jo 

- [ ds f dpe-'^*-p'^^°Ve-'^P-'^^PVe-''^° . 

Jo Js 
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